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We investigate the relation between various statistical ensembles of finite systems. If ensembles 
differ at the level of fluctuations of the order parameter, we show that the equations of states can 
present major differences. A sufficient condition for this inequivalence to survive at the thermody- 
namical limit is worked out. If energy consists in a kinetic and a potential part, the microcanonical 
ensemble does not converge towards the canonical ensemble when the partial heat capacities per 
particle fulfill the relation c^ 1 + c~ 1 < 0. 
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In most textbooks the equivalence between the differ- 
ent statistical ensembles is demonstrated at the thermo- 
dynamical limit through the Van Hove theorem 0]. In- 
deed ensembles differ at the level of fluctuations which 
are generally believed to induce small corrections in fi- 
nite systems and to become negligeable at the limit of 
infinite systems. 

In this paper we will show that this might not be al- 
ways the case. For finite systems, two ensembles which 
put different constraints on the fluctuations of the or- 
der parameter lead to qualitatively different equations of 
states close to a first order phase transition. As an ex- 
ample the microcanonical heat capacity may diverge to 
become negative while the canonical one remains always 
positive and finite [^|, [|. Such inequivalences may sur- 
vive at the thermodynamical limit for systems involving 
long range forces || . Looking at the general properties 
of the order parameter distribution a sufficient condition 
for this behavior to show up can be explitely worked out. 

Let us first concentrate on finite systems. For simplic- 
ity we will consider the microcanonical and the canonical 
ensemble characterized by the energy E and the temper- 
ature /3 _1 respectively, but our discussion is valid for any 
couple of conjugated extensive and intensive variables. 

The microcanonical ensemble is characterized by the 
level density W(E) and the entropy S — logW .The 
caloric curve is then T _1 = 8eS. The canonical par- 
tition sum is the Laplace transform of W: Zg 
j dEW exp(—/3E). In this article we will assume that 
the partition sum converges; this is not always the case 
H and indeed the impossibility to normalize the distri- 
bution Wexp(— (3E) is already a known case of ensemble 
inequivalence. 

In finite systems, the canonical ensemble differs from 
the microcanonical one since it does not correspond to a 
unique energy but to a distribution Pg(E) = cxp(S(E) — 
f3E — log Zfj). If Pg has a single maximum the average 
energy (E) B = — dg hi Zg can also be computed using 
a saddle point approximation around the most probable 
energy Eg 
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with gp(x) = C0+C3X 3 + C4X 4 + . . . . If Pp is symmetric, 
(E)g = Ep. The definition of saddle implies 
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meaning that the microcanonical caloric curve T(E) ex- 
actly coincides with the canonical one f3~ 1 ((E)). However 
in a finite system the distribution may be not symmetric 
so that the two curves can be shifted : (25) a = Eg + Sp, 

where Sg = J dx xexp(— x 2 /2C)gg(x) = 3c3\/27rC 5 + ... 
with gg the series of the odd terms of gg. However, the 
shift S is in most cases small so that when Pg has a unique 
maximum the ensembles are almost equivalent even for 
a finite system. 

A more interesting situation occurs in first order phase 
transitions where Pg has a characteristic bimodal shape 

|Q, ||, ^| with two maxima Ep, E^ that can be asso- 
ciated with the two phases and a minimum £ . These 
three solutions of Eq.(||) imply a backbending for the 
microcanonical caloric curve. A single saddle point ap- 
proximation is not valid in this case; however it is always 
possible to write Pg = m ( 3 1] ' P^ + m { g 2) ' Pf ] with P { g l) 
mono-modal normalized probability distribution peaked 

at Eg ^ . The canonical mean energy is then the weighted 
average of the two energies 
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equality holding for symmetric distributions Pg ■ As be- 
fore correcting terms depending on the skewness C3 can 
be easily derived. 

Since only one mean energy is associated with a 
given temperature the canonical caloric curve is 
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monotonous, meaning that in the first order phase tran- 
sition region the two ensembles are not equivalent 

If instead of looking at the average (E) ^ we look at 

the most probable energy Ep , this (unusual) canonical 
caloric curve is identical to the microcanonical one (see 
eq. (||)) up to the transition temperature P^ 1 for which 
the two components of Pp (E) have the same height. 
At this point the most probable energy jumps from the 
low to the high energy branch of the microcanonical 
caloric curve. The canonical curve is still monotonic 
and presents a plateau at /3 t _1 which is equivalent to the 
Maxwell construction since 
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The question arises whether this violation of ensem- 
ble equivalence survives towards the thermodynamical 
limit. This limit can be expressed as the fact that the 
thermodynamical potentials per particle converge when 
the number of particles N goes to infinity : /jv,/3 = 
/T 1 logZp/N -► fp and s N (e) = S{E)/N -> s(e) where 
e = E/N. Let us also introduce the reduced proba- 
bility pN,fi{e) = (Pp(N,E)) 1 ^ N which then converges 
towards an asymptotic distribution pn,/3 (e) — > pp (e) 
where pp (e) = exp (s(e) — (3e + fp). Since Pp(N 7 E) w 

(pp (e)) N one can see that when pp (e) is normal the rel- 
ative energy fluctuation in Pp (N, E) is suppressed by a 
factor 1/vW. At the thermodynamical limit P/3 reduces 
to a (^-function and the ensemble equivalence is recovered. 

The situation is more complicated in the case of a first 
order phase transition, i.e. for a bimodal pN,p{e) ■ As 
before, let us introduce f3Jf\ the temperature for which 
the two maxima of pn,/3 (e) have the same height. For a 
first order phase transition converges to a fixed point 

/3 t _1 as well as the two maximum energies e 



N,j3 



For all temperature lower (higher) than j3^~ only the low 
(high) energy peak will survive at the thermodynamical 
limit since the difference of the two maximum probabili- 
ties will be raised to the power N. Therefore, below ei 

and above the canonical caloric curve coincides with 
the microcanonical one in the thermodynamical limit. In 
the canonical ensemble the temperature /3 t _1 corresponds 
to a discontinuity in the state energy irrespectively of the 
behavior of the entropy between and ejp . 

The microcanonical caloric curve in the phase tran- 
sition region may either converge towards the Maxwell 
construction or keep a backbending behavior, since a 
negative heat capacity system can be thermodynamically 
stable even in the thermodynamical limit if it is isolated 
p0| . This point has been recently made in somewhat 
different words by Leyvraz and Ruffo|ll[. Examples of a 
backbending behavior at the thermodynamical limit have 
been reported for a model many-body interaction taken 
as a functional of the hypergeometric radius in ref.Q 
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FIG. 1: Temperature as a function of the potential energy 
i?2(full lines) and of the kinetic energy E — E2 (dashed lines) 
for two model equation of states of classical systems showing 
a first order phase transition. Symbols: temperatures ex- 
tracted from the most probable kinetic energy thermometer 
from eq.(5). 



and for the long range Ising model ||. This can be un- 
derstood as a general effect of long range interactions 
for which the topological anomaly leading to the convex 
intruder in the entropy is not cured by increasing the 
number of particles |l2) . Conversely, for short range 
interactions j| the backbending is a surface effect which 
should disappear at the thermodynamical limit. This is 
the case for the microcanonical model of fragmentation 
of atomic clusters [Oi and for the lattice gas model with 
fluctuating volumc|[14| . The interphase surface entropy 
goes to zero as N — » oo in these models leading to a linear 
increase of the entropy in agreement with the canonical 
predictions. From these examples, we can conclude that 
in the coexistence region the microcanonical equation of 
states may remain different from the canonical one even 
at the thermodynamical limit if the involved phenomena 
are not reduced to short range effects. 

An especially interesting situation occurs for hamilto- 
nians containing a kinetic energy contribution: if the ki- 
netic heat capacity is large enough we will now show that 
the microcanonical curve presents at the thermodynam- 
ical limit a temperature jump in complete disagreement 
with the canonical ensemble. 

Let us consider a finite system for which the hamil- 
tonian can be separated into two components E = 
Ei+E 2 , that are statistically independent (W(Ei,E 2 ) = 
Wi(Ei)W2(E2)) and such that the associated degrees of 
freedom scale in the same way with the number of parti- 
cles; we will also consider the case where Si = log W\ has 
no anomaly while S2 = log W% presents a convex intruder 
which is preserved at the thermodynamical limit. Typ- 
ical examples of E\ are given by the kinetic energy for 
a classical system with velocity independent interactions 
or other similar one body operators j^j. 

The probability to get a partial energy Eiwhen the 
total energy is E is given by 



P E (Pi) = exp (Si (P x ) + S 2 (E - E{) - S (E)) (5) 
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FIG. 2: Symbols: total caloric curves obtained with the po- 
tential and kinetic equations of state of figure 1. Lines: ther- 
modynamical limit from a (double) saddle point approxima- 
tion of the partial energy microcanonical distributions. 
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FIG. 3: Canonical event distributions in the potential versus 
kinetic energy plane (upper part) and total versus kinetic en- 
ergy plane (lower part) at the transition temperature for the 
two model equations of state of figure 1. 



The extremum of Pg is obtained for the partition- 
ing of the total energy E between the kinetic and poten- 
tial components that equalizes the two partial tempera- 
tures T = 8e 1 Si(Ei) = 8e 2 S2(E — Ei). If £ais unique, 
Pe (Ei) is mono-modal and we can use a saddle point 
approximation around this solution to compute the en- 
tropy S (E) = dE 1 exp (Si {E x ) + S 2 (E- E,)) . At 
the lowest order this leads to the microcanonical tem- 
perature of the global system 8eS(E) = T meaning 
that the most probable partial energy Exacts as a mi- 
crocanonical thermometer. If E%is always unique, the 
kinetic thermometer in the backbending region will fol- 
low the whole decrease of temperature as the total energy 
increases. Therefore, the total caloric curve will present 
the same anomaly as the potential one. 

If conversely the partial energy distribution is double 
humped Jl5| , then the equality of the partial temper- 
atures admits three solutions one of them E° being 
a minimum. At this point the partial heat capacities 
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E l ) fulfill the relation 

cr 1 



C~ L < 



(G) 



This happens when the potential heat capacity is neg- 
ative and the kinetic energy has a sufficient number of 
degrees of freedom (C\ > — C2) to act as an approximate 
heat bath: the partial energy distribution Pe (E±) in the 
microcanonical ensemble is then bimodal as the total en- 
ergy distribution Pp (E) in the canonical ensemble. In 
this case the microcanonical temperature is given by a 
weighted average of the two estimations from the two 
maxima of the kinetic energy distribution 

T = d B S(E) = f- / (7) 



where T W =T 1 (^ ) ) 



are the kinetic temperatures calcu- 



lated at the two maxima, P y "' = Pe{E 1 ) are the prob- 
abilities of the two peaks and a^ 1 ' their widths. At the 



thermodynamical limit eq.(^|) reads c^ 1 + c^ 1 < 0, with 
c = limAr^oo C/N . If this condition is fulfilled the prob- 
ability distribution Pp (E) presents two maxima for all fi- 
nite sizes and only the highest peak survives at N = 00. 

Let E t be the energy at which P Et {E {1) ) = P Et (E {2) ). 
Because of the zero principle eq.(0) at the thermodynam- 
ical limit the caloric curve will follow the high (low) en- 
ergy maximum of Pe (E\) for all energies below (above) 
E t ; there will be a temperature jump at the transition 
energy E t . 

Let us illustrate the above results with two examples 
for a classical gas of interacting particles. For the ki- 
netic energy contribution we have Si(E) = c\ \n(E/N) N 
with a constant kinetic heat capacity per particle c\ = 
3/2. For the potential part we will take two polynomial 
parametrizations of the interaction caloric curve present- 
ing a back bending which are displayed in figure 1. If 
the decrease of the partial temperature T2(^2) is steeper 
than —2/3 (left part) (J] eq.(^) is verified and the ki- 
netic caloric curve T\ (E — E\) (dashed line) crosses the 
potential one T^iE-i) (full line) in three different points 
for all values of the total energy lying inside the coexis- 
tence region. The resulting caloric curve for the whole 
system is shown in figure 2 (symbols) together with the 
thermodynamical limit (lines) evaluated from the double 
saddle point approximation (uj). In this case one observes 
a temperature jump at the transition energy while if the 
temperature decrease is smoother (right part of figures 1 
and 2) the shape of the interaction caloric curve is pre- 
served at the thermodynamical limit. 

The occurrence of a temperature jump in the thermo- 
dynamical limit is easily spotted by looking at the bidi- 
mensional canonical event distribution Pp(Ei, E2) in the 
partial energies plane. This density of states is just the 
product of the independent kinetic and potential canon- 
ical probabilities as shown in the upper part of figure 3 
for the two model equation of states of figure 1 at the 
transition temperature /? = /?(. In order to discuss the 
microcanonical ensemble one has to introduce the total 
energy E = E\ + E2 ■ Thus we can look at the canonical 
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FIG. 4: Microcanonical kinetic energy distributions inside 
the coexistence region for the two model equations of state of 
figure 1. 

distribution as a function of E and E\ 

Pp(E,Ex) ocexp5i(^i)exp5 2 (S-Si)cxp(-/3 J B) (8) 

which is shown in the lower part of figure 3. The de- 
formation of the event distribution induced by the mi- 
crocanonical constraint does not cause a topological dif- 
ference between our two model cases; this explains why 
both converge to the Maxwell construction for N — ► oo 
in the canonical ensemble. If we now study the micro- 
canonical ensemble we have to look at constant energy 
cuts of Pp(E,E\) leading to the microcanonical distri- 
bution Pe{E\) within a normalization constant. If the 
anomaly in the potential equation of state is sufficiently 
important, the distortion of events is such that one can 
still see the two phases coexist even after a sorting in en- 
ergy. Figure 4 shows two cuts of the lower part of figure 3 
at an energy close to the transition energy. In the system 
which does not present a temperature jump the partial 



energy distribution is normal, while the two peaks of the 
system violating ensemble equivalence can be interpreted 
as the precursors of phases || . This can also be seen 
from the most probable kinetic temperature (symbols in 
the left part of figure 1) which makes a sudden jump at 
the transition energy. 

In conclusion, in this paper we have analyzed the dis- 
tribution of observable quantities related to the order pa- 
rameter in finite systems undergoing a first order phase 
transition. This phenomenon is uniquely signed by a 
bimodality of the probability distribution in the inten- 
sive ensemble where the Lagrange parameter associated 
to the considered observable is fixed. In such a phys- 
ical situation the different statistical ensembles are not 
in general equivalent even at the thermodynamical limit. 
In finite systems, the order parameters (e.g. energy in 
the considered examples) have an average value which 
varies smoothly while the most probable value makes a 
jump as a function of the associated Lagrange multiplier 
(e.g. temperature) in the intensive (e.g. canonical) en- 
semble. In the corresponding extensive (e.g. microcanon- 
ical) ensemble the corresponding equation of states (e.g. 
the caloric curve) presents a back bending. In infinite 
systems, this inequivalence between statistical ensembles 
may remain. We have shown that a generic behavior of 
the extensive ensemble can be a discontinuity in the as- 
sociated intensive variable at a given value of the fixed 
extensive variable. The condition for this ensemble in- 
equivalence can be explicitly worked out in a wide range 
of physical systems. In particular, microcanonical caloric 
curves present a temperature jump at the thermodynam- 
ical limit if the negative heat capacity is sufficiently small 
in absolute value for the kinetic energy to play the role 
of a heat bath. 
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